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Abstract 

Previous results on fermion chirality-flipping four-point functions are ex- 
tended to SU (N) gauge theories. The problem is purely non-perturbative, and 
it is approached by truncating the Schwinger-Dyson hierarchy. The large-./V 
limit also simplifies the problem substantially. The resulting equation is solved 
numerically by relaxation techniques and an estimate of the critical coupling 
and momentum behavior is obtained. We also comment on the behavior of 
chirality-flipping 2n-point functions for general n. 
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I Introduction 



In this paper we study the dynamical generation of momentum dependent fermion 
four-point functions within the framework of non-abelian gauge theories. These will 
be purely non-perturbative quantities associated with breakdown of chiral symme- 
tries, and would for example imply the existence of the corresponding four-fermion 
condensates. Our analysis is based on the Schwinger-Dyson (SD) formalism. The 
resulting equations are analytically intractable and thus we seek numerical solutions. 

We have two main goals in this study. One is to estimate the critical value 
of the gauge coupling necessary for the formation of the four-point functions. We 
would like to study this in a case where four-point functions develop on scales higher 
than the scale of mass formation. This would make it consistent to treat the four- 
point function problem in isolation, independently of the mass generation problem. 
This is not the case for QCD, but even in that case it is useful to at least consider 
the existence of the effects we study here. 

In extensions of the standard model the hierarchical symmetry breaking pat- 
tern we envisage can be a natural consequence of the fact that a dynamical four-point 
function may break fewer continuous symmetries than a dynamical mass. In the 
case that some of these symmetries are gauged, then these additional interactions 
could resist the formation of mass. If so, dynamical masses may form on a scale 
lower than the scale of the four-point functions, or not form at all. With these 
possibilities in mind, we shall simply ignore fermion mass in the present work. 

A multiflavor case will be discussed elsewhere [Q] in which it is possible to add 
abelian gauge interactions such as to make the theory chiral, so that any mass will 
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break some gauge symmetry. In that theory some four-point functions do not break 
any gauge symmetry; this provides a dynamical reason for why four-point functions 
form and masses do not. There are also four-point functions for which the effects 
of the extra C/(l)'s essentially cancel out, leaving only the attractive non-abelian 
interactions. The study of such four-point functions would be very similar to the 
study of four-point functions in our simplified single-flavor case, where we assume 
the mass vanishes. 

Our other goal is to extract the momentum dependence of the dynamical 
fermion four-point function. For example if such an object is to play a role in 
generating a fermion or a gauge-boson mass, then two or four lines of the four- 
point function may be closed off into a loop or into another four-point function. 
It proves then interesting to study the details of the momentum dependence, such 
as the relative size of the four-point function when different pairs of momenta are 
large. The factorization hypothesis, according to which the four-point function can 
be treated as a product of two two-point functions, is of no use here since we neglect 
two-fermion condensates. 

In the present study we use a one gauge-boson exchange approximation. 
Since the gauge boson can attach to any pair of the four legs, the SD equation sums 
up a much more complicated set of diagrams than the set of ladder graphs. The 
present work constitutes a clear progress with respect to our previous study ||, 
since, apart from considering a general non-abelian group, it includes a treatment 
of non-linearities, it does not neglect terms proportional to external momenta while 
at the same time exploring the full available momentum space, and in a certain 
limit allows us to comment on the behavior of certain 2n-point functions. 
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In the next section we consider the SD equation for the fermion chirality- 
fiipping four-point function in a non-abelian theory, and then compute the form it 
takes in the large- N limit of an SU(N) theory in section 3. In section 4 we discuss 
our numerical treatment and results, and in section 5 we generalize our results to 
2n-point functions. 

II The equations in the non-abelian case 

We focus our attention on fermion operators which are purely chirality changing 
of the form ^l^^r^lX'^r + h.c, where tpL = ^ 2 75 ^ an d 4>r = ^ 1H 2 75 'V - We 
constrain our study to just one fermion flavor in a representation of a gauged simple 
Lie group. The four independent operators which have this property and respect 
parity areQ 

-(^^ + ^7 5 V'V'7 5 V') = 4>Lif}R$LipR+h.c. 

^a^^a^tp = ^l^^r^l^^r + h.c. 

-(ViA a ^A> + ^A a 7 5 #A a 7 5 ^) = V^A a VwALA a V^ + h.c. 

^AV'^AV'V = ij L Wip R ij L Wip R + h.c., (1) 

where A a are the generators of the non-abelian Lie algebra. The four-point functions 
corresponding to these four operators receive exclusively non-perturbative contri- 
butions. From now on we will denote all operators of the form tpTip^pT'ip, where 
r,r' are matrices with possibly non-trivial spinor and color structure, by T (g> T' . 

The chiral U(1)a is anomalous, but its discrete subgroup e i ( n / 2 ) 7r 75 w here n 

is an integer is still present and would have to be broken before fermion mass could 

x We choose to work in Euclidean space, so there is no difference between upper and lower 
Lorentz indices. 
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form. As mentioned in the introduction, other effects in a more realistic theory 
could resist the breaking of such a symmetry. Here we shall simply assume that 
fermion masses are smaller than the typical scale of our problem, and can thus be 
ignored. 

Note that we have omitted chirality flipping operators with derivatives, hav- 
ing structures like d^dj a pv <8> (1 or d^dJa pT ) and d^a pv <g> dj(j pu , where each of the 
derivatives with indices i,j, k, I = 1, 4 acts on one of the four fermion fields. We 
will present arguments supporting their omission later. 

The SD formalism relevant to these four-point functions was discussed in 
Ref. ||] and led to an equation shown diagrammatically in Fig. |l|. It results from 
a simple truncation of the SD hierarchy achieved by approximating the five-point 
function involving four fermions and a gluon by the sum of tree graphs involving the 
four-point function. (The equation has been symmetrized to include diagrams with 
gluons connecting all possible pairs of fermions, so the right-hand side is multiplied 
by a factor of 1/2.) The analogous procedure used in the SD equation for the 
two-point function yields the popular ladder approximation. In our case the set 
of diagrams being resummed has a more complicated structure than sets of ladder 
graphs. In both cases gauge self-interactions no longer appear in the SD equation. 
In the two-point case it is then true that there is little difference between abelian 
and non-abelian interactions. On the contrary in the four-point case we show that 
there is still a big difference; for non-abelian interactions the space of four point 
functions is enlarged and the mixing between the four-point functions is completely 
different. 

We are considering a four-point function associated with the Green function 
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(O|T{ , q ,V'/3V'pV't}|O). In momentum space, the four-point function involving scalar 
functions of momentum only and receiving exclusively non-perturbative contribu- 
tions is 

+O c s °l°£(\ a I a p ® A% T + X a 7 5 a/3 ® A a 7p 5 T ) + Of° r \ a o% ) ® AV^ (2) 

where the four scalar functions Os+p, Ot, O^tp and 0^ lor depend on 6 variables, 
which are all the independent and Lorentz-invariant combinations of the external 
4-momenta p\, ...,^4, corresponding to the fermions with spinor indices a,(3,p,r 
respectively. We wish to develop the SD equations for these scalar functions. 

Similarly to our previous work Q, we begin by defining the following func- 
tional operators T l : 



8vr 3 J (pi - k) 2 (p 2 + k) 

r»[Jt] s £ ( 
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8vr 3 7 (p 3 + k) 2 (p4- k) 
T c lK] s JL[, k . * 
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87r 3 7 (Pi-A;) 2 (p4-A:) 
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8vr 3 7 (p 2 + A;) 2 (p 3 + ^) 
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8vr 3 7 (Pi-A;) 2 (p3-A:) 2 

rF[K] = dAk {P2+k f(p 4 +kr (3) 

where X is a function of the loop and external momenta, with a possibly non-trivial 
spinor structure, the letters A, F correspond to the diagrams shown in Fig.^J with 
the gauge boson having four-momentum k, and a is the momentum-independent 
coupling. We work in the Landau gauge which is popular in studies of SD equations; 
the gauge boson propagator reads ^p- = p- ( — ^jr- ) • The choice of this gauge 



will be further justified later, when we consider the large- A limit of an SU(N) gauge 
theory. 



By combining the results of the abelian case Q with the study on the color 
structure given in Appendix A, and considering fermions in the fundamental repre- 
sentation of SU(N), we have 

Os+p = 3{N '~ 1} (f* + r*) [O s+ p] +6(fc + fD-fE-fF) pf"] 



o T = - ^ 1 (r^ + T b )[O t ] + 2{T C + T D + Y E + T F )[Of lor ] + 



I( F C + T D _ t e _ r F )[Og$Z] 



0g& = -^(r A + rg)[ogjg] + 3{N 2N2 1} (ro + r*> - v E - t*)[Ot\ + 



/ /V 2 -I- 9 1 \ 



Ofor = ^ (r A + rB)+ ^_^^ (r C + r D ) _^ (rS + r i, ) jj [0 co^ ] + 



A 2 — 1 -TT —FT -tt, r „ , A 2 - 1 



2 y , (re + + + v E )[o T ] + ^^(r c + - - r E )[o s+P } + 

1 / /V 2 -I- 9 1 \ 

2 ( ^ + + N {TE + TF) ) (4) 

where T % [0] = r l (0 l ) and O l are the form factors with loop-dependent arguments 
corresponding to diagrams i = A, ...,F 0. 
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We are faced with a system of four coupled 4-dimensional integral equa- 
tions involving functions of 6 variables. As it stands, the problem is analytically 
intractable, and even a numerical solution proves to be beyond our present means. 
In the next section we present the simplification that the large- N limit of SU(N) 
provides. 

Ill The large- TV limit 

By taking the large- N limit, (||) reduces to [] 



( O s+ P \ 



sncolor 
U S+P 



N 



( |(rA + ffl) 
o 
o 



V 









.l(fX + fB) 















-3(rc + fi>) 
_(fc + f^) j 



( [Os+p] \ 
Pt] 



U S+P 



Qcolor 



The form factors Os+p and Op decouple from the others and from each other. The 
Os+p entry is positive, which means that there is an attractive interaction in this 
channel necessary for the formation of a non-zero four-point function. In contrast 
Op is not expected to be non-zero. It is also apparent that O c s ° + °p can be non-zero 
only if Oj? lor is somehow non-zero, but that seems unlikely because of the negative 
entry coming from Of lor itself. Even if the O c s ° + % T functions were somehow non- 
zero, we would expect for them a critical coupling larger than the one required for 

2 Note that we apply the large- JV limit only within the context of the truncated SD equation, 
and we do not claim that the large- N limit in any way justifies the original truncation of the SD 
hierarchy. 



(5) 
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Os+p- 

The large- N limit also allows us to argue for the omission of contributions 
coming from operators proportional to external momenta. With regards to the 
operators with structure p^a^ ®pja uX and pfpjCr^ 11 ^>p^pJcr XT , we can see that the 
contributions they receive from the dominant diagrams having the same functions 
in their vertices (diagrams A and B for 1 (g> 1 color structures, diagrams C and D for 
A a <g> A a color structures) are negative, so we do not expect the corresponding four- 
point functions to be non-zero. This happens because, in the large-iV limit, tensor 
insertions give negative contributions. Even if other terms made them non-zero, 
they would require a larger gauge coupling than the one needed for Os+p- 

Interesting also are the operators with structure 0\ = p^p^a^ <g> 1 and 
O2 = 1 <8> pfpjCr^ (and trivial color structure), which can obviously mix with 
the form factor Os+p- If we take the operator 0\ for instance, we see that it 
receives a negative contribution from diagram A and a positive one from diagram 
B. Therefore, it might not develop a non-zero value, and if it did, it would require 
a larger value for the gauge coupling than the critical value we are seeking. These 
operators will also be neglected. 

We now write down the integral equation for Os+p- We go to a reference 
frame where p\ = —p~2 , where the form of the kernel of the integral equation indicates 
that Os+p is a function of five variables, i.e. 

x f Os+p(k°,\k\,pl \p 4 \, g°)D^ | O s+ pip i,\p'i\,k ,\k\,Q )Dl 
\ (k-pi) 2 (k-p A ) 2 

where the two terms on the right-hand side come from diagrams A and B respec- 
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tively, q = Pi + P2, D^ v = — - — P (j~_ p .y2 > an d the two fermions in the loop 
have four-momenta k and k — q respectively. 

We find that a solution to the above integral equation takes the factorized 
form Os+p ~ B(pi, \p{\, q )B(p®, \pl\,q ). This just reflects the fact that the loop 
integral in diagram A does not depend on and \p^\ for a given q°, and similarly 
for diagram B. By momentarily setting p\ = p\ = p° and \p~[\ = \p~l\ = \p\ we obtain 
an equation for B. 



We have separated q from the other two variables with a semi-colon since it is an 
argument of B that is not affected by loop momenta, and it can thus be treated as 
if it were a parameter in the kernel. 

Note that for q° = this reduces to the linearized SD equation for the fermion 
two-point function. For non-zero q°, the problem reduces to the one of a three-point 
function, which has been studied before ||. After some Dirac algebra, and omitting 
terms with structure p^p^a^ v 1 according to our previous discussion, the equation 



In the following we will change our variables from p°, k°, \p\, and |Ar|, to p = 



y/pO 2 _|_ _ 2 _|_ |£|2 ; ^ — arc tan (\p\/p°) and <f> = arctan (\k\/k°). With- 

out loss of generality we will write B(p°, \p\;q°) = R(q°)B(p,(j);q ) with R(0) = 1 
and B{0,0;q°) = 1. 




becomes 
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At this stage our equation is linear, and the function R(q°) remains unde- 
termined. Clearly some essential physics, containing the non-linearities, has been 
omitted by the original truncation. In the two-point function case it is known that 
the main effect of the non-linearities can be modeled by introducing an infrared 
cutoff, which is determined in a self-consistent manner by identifying it with the 
dynamical mass evaluated at the infrared cutoff. We will follow a similar procedure 
in the four-point function case by identifying the infrared cutoff with AR(q°), where 
A is the infrared cutoff when q° = 0, and then determining R(q°) in a self-consistent 
manner. 

Non-linear effects will enter through diagrams involving an odd number of 
four-point functions, an example of which is provided by the diagram in Fig.||. For 
large gluon momentum this diagram is much more suppressed than the diagrams 
we have been considering; the diagram only becomes important for small gluon 
momentum. This diagram thus provides an example of how non-linear effects could 
introduce a natural infrared cutoff into our problem. Unfortunately, it is not easy 
to resum such diagrams, as is possible in the two-point function case. 

Our result for the 4-point function will take the form 

O s+P = R^) 2 B{pi,<h.^)BMi\q°), (9) 

where r is a dimensionless constant. It is interesting that r may have to be signifi- 
cantly larger than unity, in order to compensate for the loop-phase-space suppression 
factors appearing in diagrams such as Fig.^. 

In place of an explicit infrared cutoff we choose to multiply the kernel of our 
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equation by the function k /(k + A R (q )). The final equation becomes 

„, : 3iVa /■ ~ kB(k,4>;q°)k sin 4> 
pB(p,(f);q ) = —r- \ dhd^-T— — — =r X 



471-2 ^ ' - kcos<t>) 2 + k 2 sm 2 4>^j 2sin0 



/ ( A: sin (p + p sin </>) 2 + (fc cos (p — p cos (/>) 2 \ 
y (fc sin </> — p sin <?!>) 2 + (A; cos — p cos 0) 2 / 



/ q°k cos (f>\ k 2 

y l~ 2 ) x k 2 + K 2 R 2 { q ^y { ' 

where f 0, f vr. We solve for the function pi? in order to increase the 
accuracy of the numerical solution that follows, since we expect B to decrease with 
increasing p. We note that our crude representation of the infrared effects does not 
affect the functional form of B when p is large. 

With regards to the gauge dependence of our results we note that the use 
of the Landau gauge, along with a bare massless fermion propagator and a bare 
gluon-fermion vertex, is consistent with the Ward-Takahashi identity at one-loop. 
Moreover, in the large- A 7 limit only diagrams A and B are important and the situ- 
ation is then expected to be similar to that of the two-point function. In particular 
in that case the use of the Landau gauge in the ladder approximation yields results 
resembling those found in gauge-invariant treatments j|. 

IV Numerical results 

The form of the integral equation allows us to use the same discretization lattice 
for the arguments of the function B inside and outside the integral. The angles 
are discretized according to 4>{i), 4>{i) = ' ^ = ^' •••' n - We- avoid the endpoints 
where the angles are equal to or ir since our integral has an integrable singularity 
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there. 

The momenta are discretized according to log w p(i),log 10 k(i) = log 10 (A/.R + 
^j-log 10 (A[/v//Aj#), where h-iR^uv are the infrared and ultraviolet cut-offs. Ajjv 
does not have to be much larger than the energy scale where our function is negli- 
gibly small. On the other hand, Am has to be smaller than A which provides an 
effective infrared cutoff to our theory. We choose to work with the ratios A/ Am = 10 
and Au\r/A = 100. We choose the number of points (n) in each dimension to be 
equal to 40. Our results change little when this number is increased. 

The integral equation is solved via a relaxation method. We first insert 
an initial configuration for our function, and then iterate the equation until it is 
satisfied to a reasonable accuracy. We start by setting q° = and determine the 
critical coupling necessary for a solution. We then keep this coupling fixed, and 
for different values of q° we compute B(p,(/);q ) and R(q°). Our critical coupling 
is Na c = 2.7 ± 0.3 which, from our previous discussion, is the same as the one for 
chiral symmetry breaking via the two-point function with the same cut-offs. In the 
large- N limit and for infinite Auv this coupling is given by Na c = ^ ~ 2.1. 

The function B(p,4>;q°) has only a weak dependence on q°, and in Fig. || 
we plot j^B(p, (j>; 0). We see that j^B(p, cj>; 0) exhibits the cos (7 log (p/A)) behavior 
(with 7 a coupling-dependent constant) that is well known in two-point function 
studies. B(p,(ft;0) is independent of <fi as expected, and even for non-zero q° the 
dependence on <f> remains very weak. 

The main q° dependence enters in the function R(q°) which we plot in Fig.|j. 
R(q°) falls rapidly with increasing q°; we find that q° cannot exceed A by much 
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(ps 1.16A) before the whole solution collapses to zero. 

V 2n-point Green functions 

The form that the integral equation takes for large N suggests a possible gener- 
alization of the above results. In particular, the diagrams that dominate are the 
ones in which a gluon is attached between two fermion fields with spinor indices 
contracted with each other. Therefore, one might expect only condensates of the 
form < ^l^r-^l^r + h.c. >, where the ellipsis stands for 2n — 4 fermion fields 
paired with each other. They should form at a critical value of the gauge coupling 
similar to that required by the two- and four-point functions. 

Moreover, we expect the corresponding scalar 2n-point function Og +P to 
factorize as Og +P ~ R n B\...B n . Each B- L depends on the momenta of the respective 
fermion pair. In order for each of the Bi to satisfy the same final integral equation as 
before, however, we would have to assume that the three-momenta of the fermions 
in each pair are equal and opposite, so that ql = 0. We expect, though, that our 
previous functions B and R should describe the qualitative momentum behavior of 
the 2n-point functions if we are able to go to an optimal reference frame where the 
various ql are reasonably small. 

VI Conclusions 

In this work we considered chirality-changing fermion four-point functions in a non- 
abelian gauge theory which receive non-perturbative contributions exclusively. We 
have tackled the problem by writing down an equation derived from a truncation 
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of the infinite SD hierarchy. The solution to this equation is expected to illustrate 
some generic properties of four-point functions. The large-iV limit of SU(N) ren- 
ders the four-point problem similar to that of the three-point function. We find 
that a particular four-point function is much more likely to form than other four- 
point functions, and that its critical coupling is similar to the one required for a 
two-fermion condensate. The large- N limit also allows us to discuss the critical and 
momentum behavior of fermion chirality-flipping 2n-point functions for n higher 
than two. 
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Appendix A 

The purpose of this appendix is to investigate the effects that the gluons in dia- 
grams A, ...,F have on the color structure of the fermion vertices. We consider the 
case of fermions in a general representation of a compact simple Lie group. The 
group is then defined by the commutation relations [A a ,A 6 ] = t abc X c , where A a are 
the (traceless) generators of the corresponding algebra. 

The Dynkin index T f is defined by the relation Tr(X a X b ) = T f 5 ab . By N f = 
Trl we denote the dimension of the fermion representation, and by Cf and C g the 
quadratic Casimirs of the fermion and adjoint representations. These are defined by 
A a A a = Cfl and t abc t abc ' = C g 5 cc ' . For fermions in the fundamental representation 
of SU(N), N f = C 9 = N, T f = I, and C f = 
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There are two color structures in the problem. One is l(g>l, which is associ- 
ated with the functions Os+p,t, and the other is A a <g> A a , which is associated with 
the functions Of^ T . 

I. VERTEX COLOR STRUCTURE: 1® 1 

i. Diagrams A, B 

We have A a A a <8> 1 = C/l ® 1. Diagrams multiplied by C/, no change in color 
structure. 

ii. Diagrams C,D,E,F 

Color structure becomes A a <8> A a . 

II. VERTEX COLOR STRUCTURE: A a <g> A a 

i. Diagrams A, B 

X b X a X b ®\ a = (C f - C g /2)\ a ® X a . Diagrams multiplied by C f - C g /2, no change 
in color structure. 

ii. Diagrams E, F 

Color structure becomes Oef = A a A fe <8) A a A 6 . This has to be reduced down to the 
two initial color structures. 

We first expand the matrix A a A fe with respect to the basis of the unity matrix 
1 and the group generators A c . 

X a x b = f a bl + h abc x c (n) 

By taking the trace of both sides, we find that f ab = j^j5 ah . Moreover, considering 
the case a = b proves that h aac = 0. 
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The above considerations allow us to write 

Oef = |fc / l<g>l + h abc h abc '\ c ® A c ' (12) 

We now have to calculate the quantity h abc h abc ' . By virtue of (|TT|), we con- 
sider the tensor 

h abc h abc> x c x c> = ( _ Ifs^A f ^6 _ TLs^I] 

\ N f A ^ / 

But the unique quadratic tensor proportional to 1 is the quadratic Casimir of the 
fermion representation. This means that h abc h abc ' = h5 cc \ in which case we find 

So finally we have a new structure Oef = jf-CfX (g> 1 + /iA c ® A c , with /i given 
above. For fermions in the fundamental representation of SU(N), h = —1/N. 

iii. Diagrams C, D 

Color structure becomes Ocd = A a A b ® A 6 A a . By using the commutation relations 

C 

of the group, we find that the color structure changes to Ocd = Oef — ^-A a <8> A a . 
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Figure 1: The schematic form of the SD equation. We have labeled the four fermions 
by their spinor indices. We label the diagrams by the capital letters A, ...,F, and 
we omit the factor of 1/2 multiplying the right-hand side. 



Figure 2: An example of a diagram which introduces non-linearities in the infrared. 



Figure 3: The function %B(p, <f>\ 0) 



Figure 4: The function R(q°). 
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Figure 1: The schematic form of the SD equation. We have labeled the four fermions 
by their spinor indices. We label the diagrams by the capital letters A, ...,F, and 
we omit the factor of 1/2 multiplying the right-hand side. 




Figure 2: An example of the non-linear terms that act as an IR cut-off in our theory. 




Figure 3: The function ^BfafcO) 
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